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Introduction 

The  amount  of  water  in  the  optical  path  between  an  airborne  measurement  platform  and 
a  target  is  responsible  for  the  difference  between  the  actual  and  the  measured  target  spectral 
radiance.  In  situ  measurements  of  water  vapor  can  be  limited  by  outgassing  from  the  sensor 
and/or  its  environment.  The  objective  of  this  study  was  to  investigate  the  use  of  lidar  remote 
sensing  techniques  for  the  remote  measurement  of  water  vap)or  from  a  measurement  aircraft  such 
as  the  FISTA  NKC-135A. 

This  report  is  a  brief  study  of  the  signals  anticipated  in  measuring  the  water  vapor 
concentration  in  the  troposphere  from  a  lidar  flown  on  an  aircraft.  Two  detection  mectiamsp'S 
have  been  examined:  conventional  Raman  scattering  and  differential  absorption. 

The  basic  lidar  equation  is 

e-  A  5r/r^ 

where  P,,  =  incident  laser  signal 

€  =  system  efficiency 

e  '  =  atmospheric  transmission 
/3  =  backscatter  coefficient 
A  =  collecting  area 
5r  =  size  of  atmospheric  cell  probed 
r  =  range 

The  two  methods  of  detection  are  significantly  different  in  approach.  Raman  scattering 
focuses  on  the  backscatter  coefficient,  0,  to  obtain  specie  dependent  information.  The 
differential  absorption  technique  is  independent  of  the  scattering  mechanism  and  obtains  the 
information  from  the  extinction  term,  e  ^ 

In  both  cases  the  laser  beam  is  scattered  180  degrees  back  to  the  receiver.  Fig.  1.  shows 
the  ratio  of  /3„ro»o//3„yieigh  for  different  wavelengths  as  a  function  of  altitude.  Backscatter  from 
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aerosols  dominate  for  the  first  7  km  whereas  the  Rayleigh  scattering  is  responsible  for  scattering 
between  the  altitudes  of  7  and  15  kilometers.  These  data  have  been  fit  to  a  polynomial  and  are 
extrapolated  linearly  for  wavelengths  greater  than  1  ^m. 

Water  Vapor  in  the  Troposphere 

The  concentration  of  water  vapor  in  the  troposphere  varies  by  over  a  factor  of  10000. 
The  mixing  ratio  (ppmv)  decreases  from  7750  ppmv  at  sea  level  to  5  ppmv  at  15  km  and 
remains  essentially  constant  through  out  the  stratosphere  and  mesosphere.  The  corresponding 
number  density  of  water  molecules  is  2x10^^  m'^  and  2x1  O'*  m'^  at  sea  level  and  15  km, 
respectively.  Fig.  2.  shows  a  plot  of  the  water  vapor  density  as  a  tunction  of  altitude  over  these 
altitudes. 

Optical  System 

The  calculations  in  this  report  were  made  assuming  a  coaxial  lidar  system  in  order  to 
utilize  the  entire  near  field.  In  addition,  the  collecting  optic  is  assumed  to  have  a  7  inch 
diameter  (248  cm^  collecting  mirror.  The  optics  are  assumed  to  consist  of  the  aircraft  window, 
a  telescope,  a  filter  and  lens,  and  detector.  Below  is  a  list  of  the  efficiencies  used  for  the 
following  set  of  calculations 


Optical  Efficiencies 


Aircraft  Window  (Trans) 

0.90 

Aircraft  Window  (Rec) 

0.90 

Telescope 

0.90 

Filter 

0.50 

Lens 

0.90 

PMT  Detector 

0.10 

TOTAL 

0.033 

2 


If  an  InAs  detector  having  an  0.50  quantum  efficiency  is  used,  the  total  optical  efficiency 
becomes  0.164. 

Aircraft  Contamination 

An  added  difficulty  in  trying  to  measure  water  is  that  most  materials  will  experience 
some  outgassing  of  water  and  potentially  contaminate  the  experiment.  Simple  calculations  were 
performed  to  examine  whether  outgassed  water  from  the  aircraft  could  interfere  with  any 
measurements  being  made  from  the  plane. 

When  the  diffusion  of  water  vapor  is  considered  by  itself,  it  is  found  that  at  15  km, 
where  the  mean  free  path  of  water  vapor,  X  is  largest,  any  outgassing  of  water  from  the  aircraft 
skin  remains  very  close  to  the  skin  and  will  not  readily  contaminate  any  measurement. 

More  importantly  is  the  effect  of  a  transverse  high  altitude  wind  with  a  velocity  v^. 
Water  vapor  outgassed  from  the  aircraft  will  be  carried  along  by  the  wind  a  distance,  d  =  v^t 
where  t  =  z/v,  and  z  is  the  distance  from  the  tip  of  the  aircraft  to  the  station  of  the  sensor. 
The  distance  to  the  edge  of  the  outgassed  water  cloud  is 


using  reasonable  values 

d  =  =  2.7  10^  cm 

2.7  X  10“  cm/s 


which  is  still  much  smaller  than  typical  lidar  ranges 
Conventional  Raman  Scattering 

Raman  scattering  is  a  technique  that  identifies  target  species  by  a  unique  scattering 
mechanism.  The  incident  wavelength  is  absorbed  by  the  molecule  which  is  excited  into  a  virtual 
state.  The  scattered  wavelength  is  emitted  with  a  wavelength  that  is  different  than  the  incident 
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photon  and  the  molecule  is  left  in  some  final  state.  Raman  scattering  is  a  useful  method  for 
identifying  specific  constituents  from  a  mixture  of  gases.  The  vibrational-rotational  cross 
sections  for  this  phenomena  are  typically  several  orders  of  magnitude  smaller  than  Rayleigh 
scattering. 

Simulations  of  Raman  scattered  signals  from  water  vapor  have  been  made  for  the  Nd:  Yag 
wavelengths  532,  355,  and  266  nanometers.  The  incident  energy  per  pulse  used  for  the  lasers 
was  150,  30,  and  80  mj,  respectively.  (N.B.  The  266  nm  power  level  was  obtained  from  a  laser 
catalog;  the  532  and  355  nm  wavelengths  are  based  on  laboratory  measurements.)  The  results 
are  directly  proportional  to  the  energy  per  pulse.  Although  the  532  nm  wavelength  has  more 
incident  power  than  the  UV  wavelengths,  the  combination  of  higher  background,  dark  counts, 
and  the  X'"*  cross  section  dep)endence  make  it  a  smaller  signal  to  noise  ratio.  The  return  signals 
have  been  integrated  for  200  seconds  at  10  pps  (  2000  shots  coadded)  with  a  sampling  rate  of 
0.2  US  (30  meter  bins)  and  are  shown  plotted  in  Fig.  3a-c.  The  long  integration  time  is  required 
to  obtain  a  useful  measurement.  The  data  from  these  plots  are  also  presented  in  tabular  form 
in  Appendix  A. 

The  calculations  show  that  measurements  with  adequate  signal  strength  would  be  confined 
to  a  range  of  a  kilometer  at  low  altitudes,  hundreds  of  meters  from  the  aircraft  around  5  km,  and 
less  than  100  meters  above  10  km.  The  266  nm  signal  is  the  strongest  of  the  three  wavelengllis 
followed  by  the  355  nm  signal;  this  is  because  of  the  strong  X  '*  dependence  of  the  cross  section. 
The  signal  decreases  rapidly  for  altitudes  above  8  km  due  to  the  fall  off  of  the  mixing  ratio.  In 
all  three  cases  the  signal  above  i2  km  is  vciy  weak  to  the  i^i.craft  where  I.'r^  is 

rapidly  increasing.  Fig.  4  is  a  plot  of  range  versus  altitude  for  the  different  wavelengths.  Each 
point  on  a  curve  represents  a  return  signal  of  ~  100  photoelectrons  which  is  a  statistically 
significant  (10  f)ercent)  signal. 

Differential  Absorption 

The  differential  absorption  lidar  technique  (DIAL)  is  a  measurement  of  the  extinction 
term  and  is  independent  of  0.  Whereas  measurements  of  Raman  scattering  are  dependent  on  the 
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backscattering  of  photons,  differential  absorption  measurements  are  dependent  on  what  happens 
to  the  photons  before  and  after  backscattering.  The  extinction  coefficient  is  separated  into  two 
terms 

T  =  x(X,r)  or  +  N(r)  a{X,r) 

One  term  contains  the  effects  from  all  atmospheric  constituents  except  the  species  of  interest  and 
the  other  term  is  due  to  absorption  exclusively  by  the  target  species.  The  target  concentration 
is  determined  by  comparing  a  measurement  at  a  wavelength  sensitive  to  the  absorption  with  a 
second  measurement  very  close  in  wavelength  but  insensitive  to  absorption  by  the  target  species. 
The  two  equations  are  solved  for  the  density  which  is  given  below: 

„  _  1  [P(X,,r)  P(X,,r-<Sr)] 

2[a(X,,r)  -a(Xj,r)  ]  S  r  [P(X,,r-£r)  P(\,,r)  ] 


Choice  of  Absorption  Line 

There  are  hundreds  of  water  absorption  lines  to  choose  from.  Several  criteria  are  used 
to  determine  the  optimum  line.  The  line  should  be  chosen  to  be  in  an  spectral  region  where 
other  constituents  (e.g.  CO2)  do  not  interfere.  It  is  desirable  that  the  absorption  cross  section 
He  independent  of  temperature,  therefore,  a  ground  state  transition  is  desirable.  The  line  should 
be  chosen  such  that  the  extinction  does  not  drive  the  exponential  into  saturation  and  is  the  most 
sensitive  region  of  the  exponential. 

Laser  Spectral  Bandwidth 

The  equation  derived  for  N(r)  assumes  that  the  laser  is  monochromatic.  When  the 
spectral  bandwidth  of  the  laser  is  narrow  compared  to  the  absorption  peak  this  is  a  good 
approximation,  however,  when  the  laser  bandwidth  is  wide  (e.g.  dye  laser)  or  the  absorption 
feature  is  respectively  narrow,  the  amount  of  absorption  must  be  integrated  over  the  wavelength. 
This  requires  an  accurate  knowledge  of  the  spectral  shapes. 
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Simulations 


DIAL  returns  for  wavelengths  "otT  the  absorption  peak  were  calculated.  The  assumption 
is  made  that  an  appropriate  absorption  line  can  be  found  providing  enough  signal  can  be 
backscattered  off  the  atmosphere.  These  calculations  do  not  include  background  or  detector 
noise.  The  calculations  have  been  carried  out  for  the  altitudes  or  2  ,5  ,7  ,10  ,12  ,and  14  km 
and  have  coadded  2000  shots  in  30  meter  bins  for  direct  comparison  with  the  Raman  scattering 
data. 

The  calculations  are  made  at  1.85  ;tm,  a  region  where  several  suitable  K,0  absorption 
lines  are  found  (e.g.  X  =  1.8549  ^m,  (1=1.55x10  ’*  cm’).  A  plot  of  the  data  is  shown  in  Fig. 
5;  the  same  data  is  listed  in  a  table  in  Appendix  B.  The  return  signals  have  gcxxl  statistics  out 
to  2  -  3  km  at  all  altitudes.  Note  that  the  measurement  at  2  km  is  less  than  the  signal  at  5  km 
at  a  range  greater  than  ~  270  meters.  This  is  due  to  the  peak  in  the  backscatter  mixing  ratio 
seen  in  Fig.  1  and  the  increased  extinction  at  the  tower  altitude 

The  return  signals  at  1.85  ^m  are  significantly  larger  than  those  due  to  Raman  scattering. 
This  is  because  of  the  small  size  of  the  Raman  scattering  cross  section.  The  DIAL  method 
generates  875  photoelectrons  from  a  range  of  2  km  and  an  altitude  of  14  km.  At  a  range  of 
1000  meters,  typical  return  signals  are  in  excess  of  lOf'  photoelectrons  at  altitudes  below  7  km 
and  a  few  thousand  counts  up  to  14  km.  The  above  numbers  do  not  include  background  or 
detector  noise. 

Detection  at  1 .85  /(in  can  be  done  with  a  InAs  detector  cooled  to  LN?  temjjeratures.  The 
value  for  D’  at  this  wavelength  is  -  10”  cm-Hz''^/W  from  which  a  NF.P  of  --  3000  photons/bin 
is  obtained.  The  natural  background  was  assumed  to  be  due  to  scattering  from  a  300  K 
blackbody.  The  300  degree  blackbody  radiance  at  1.85  fim  is  9.5x10'’  Watts-cm  ’-/um.  After 
folding  in  the  system  characteristics  the  blackbody  contribution  is  negligible.  The  total 
background  is  therefore  assumed  to  be  defector  noise 

When  the  background  is  considered,  the  signal  at  the  lower  altitudes  still  maintains  gcxid 
signal  to  noise  out  to  ranges  over  a  kilometer.  In  a  worst  case,  the  range  is  restricted  to  ~400 
in  at  14  km. 

A  suitable  laser  emitting  at  this  wavelength  is  available  trenn  Schwartz.  Electro-Optics. 
The  laser  is  a  solid  state  laser  using  a  CoiMgF,  crystal  for  lasing  (Nd:YAG  pumped)  and  is 
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tunable  over  the  range  1.75  to  2.5  fxm.  The  laser  can  be  Q-switched  with  a  typical  pulse  energy 
of  20  mj  and  a  pulse  width  of  —  100  ns  (15  meters).  The  spectral  bandwidth  is  ~0.1  cm 
however  other  optics  can  be  added  to  further  reduce  the  width  and  still  maintain  the  output 
power.  This  is  because  the  output  power  is  determined  by  the  mirror  damage  threshold  and  the 
laser  can  be  pumped  harder  until  the  damage  threshold  is  reached.  A  reference  for  some  of  this 
work  is  found  in  the  IEEE  Quantum  Electronics,  Vol.  QE  21,  p  202,  1985. 

The  eyesafe  characteristics  make  the  above  laser  a  viable  candidate  for  airborne 
measurements. 

Topographical  Scatter 

One  possible  way  to  extend  the  range  of  the  measurement  is  to  replace  with  the 

backscatter  coefficient  from  a  hard  target.  This  would  significantly  increase  the  signal  to  noise 
ratio,  however,  range  information  would  be  sacrificed.  This  target  might  be  in  the  form  of  a 
large  area  retroreflector  mounted  on  the  side  of  another  aircraft.  Retroreflector  arrays  with  a 
ty^iical  reflector  dimension  of  150  ^m  are  etched  onto  0.060  inch  thick  alumina  thus  providing 
large  area  retroreflectors  which  are  light  and  compact. 

Summary  and  Recommendations 

Preliminary  sensitivity  calculations  have  been  made  for  detecting  water  vapor  in  the 
troposphere  and  show  that  it  is  feasible  to  make  remote  measurements  from  an  aircraft.  The 
concentration  of  water  vapwr  covers  a  very  large  dynamic  range  (xlO*  over  0  -  15  km). 
Conventional  Raman  scattering  at  Nd.Yag  wavelengths  and  differential  absorption  (X=^1.85^m) 
detection  mechanisms  have  been  examined. 

The  results  from  Raman  scattering  show  that  the  signals  must  be  integrated  for  a  large 
amount  of  time  (e.g.  200  seconds)  to  obtain  statistically  significant  measurements.  The 
calculations  show  that  measurements  would  be  confined  to  a  range  of  a  kilometer  at  low 
altitudes,  hundreds  of  meters  from  the  aircraft  around  5  km,  and  less  than  100  meters  above  10 
km.  Data  regarding  longer  ranges  must  be  inferred  statistically.  Scattering  at  the  shorter 
wavelengths  is  favored  because  of  the  strong  X  dependence  of  the  cross  section. 
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The  DIAL  method  is  a  more  sensitive  technique  and  therefore  a  more  promising 
technique  for  an  aircraft  lidar.  The  DIAL  technique  is  a  more  complicated  experiment  and 
requires  measurements  at  a  second  wavelength.  However,  it  is  able  to  record  measurements  on 
the  order  of  a  couple  of  kilometers  at  lower  altitudes  ( <  7  km)  and  typically  0.5  kilometer  at  the 
’  igher  altitudes.  The  use  of  a  hard  scatterer  can  significantly  enhance  the  SNR,  but  range 
information  is  sacrificed. 

Based  on  the  results  of  this  study  it  is  recommended  that  a  preliminary  engineering  design 
for  an  airborne  DIAL  lidar  be  developed.  Using  the  simulation  tools  developed  for  this  study, 
the  engineering  design  specification  will  be  optimized  and  a  preliminary  flight  test  plan 
developed. 
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BACKSCATTER  MIXING  RATIO 


Figure  1.;  Mixing  ratio,  for  various  wavelengths  (  from  Russell,  P.B.,  B.M. 

Morley,  J.M.  Livingston,  G.W.  Grams,  and  E.M.  Patterson,  "Improved  Simulation  of 
Aerosol,  Cloud,  and  Density  Measurements  by  Shuttle  Lidar",  NASA  Contractor  Report 
3473,  NASA  Langley  Research  Center,  1981) 
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Figure  3a.:  Raman  scattering  at  532  nm 
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Figure  4.:  Range  versus  altitude  for  constant  Raman  scattering  returns 


Appendix  A:  Raman  scattering  returns  at  532,  355,  a.id  266  nm 
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